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ABSTRACT: We investigate how the predictions of the Minimal Supersymmetric 
Standard Model are modified by D-term contributions to soft scalar masses, which 
arise whenever the rank of the gauge group at very high energies is greater than 
four. We give a parameterization of the most general such contributions that can 
occur when the unbroken gauge symmetry is an arbitrary subgroup of Eq, and 
show how the D-term contributions leave their imprint on physics at ordinary 
energies. We impose experimental constraints on the resulting parameter space 
and discuss some features of the resulting supersymmetric spectrum which differ 
from the predictions obtained with universal boundary conditions on scalar masses 
near the Planck scale. These include relations between squark and slepton masses; 
the behavior of sin 2 (/3 — a) (which determines the production cross-section for the 
lightest Higgs scalar boson at an e + e~ collider) and the mass of the pseudoscalar 
Higgs bosons; [the ratio T(Z — > bb)/T(Z — > hadrons)]; and mass differences 
between charginos and neutralinos. 

' Address after Sept. 1, 1995: School of Natural Sciences, Institute for Advanced Study, 
Princeton, NJ 08540 



1. Introduction 



The mass of the Higgs scalar boson in the standard model is subject to quadratically 
divergent radiative corrections of order 5m? H ~ A 2 , where A is an ultraviolet cutoff mass. 
This poses a naturalness problem for A much larger than the electroweak scale. Low- 
energy supersymmetry (SUSY) [1] presents a beautiful solution to this problem, because 
all quadratic divergences cancel order-by-order in perturbation theory. The cancellation 
still works if SUSY is softly broken by scalar mass terms, scalar trilinear couplings, and 
gaugino mass terms. These are exactly the parameters which determine the masses of 
the superpartners of standard model particles; therefore, in the minimal supersymmetric 
standard model (MSSM), the superpartners must have masses which do not greatly exceed 
the TeV range in order to maintain the stability of the electroweak scale. Moreover, we 
already know that the soft SUSY-breaking terms in the MSSM are far from arbitrary in 
their flavor and phase structure. Otherwise, large flavor changing neutral currents could be 
expected to manifest themselves in processes such as K — K mixing, b — > S7, and \i — > e^y. 
Arbitrary complex phases would give rise to a CP- violating electric dipole moment for the 
neutron in violation of experimental bounds. Thus there is strong circumstantial evidence 
in favor of some organizing principle governing the soft SUSY-breaking terms. 

Supergravity models [2] do provide just such an organizing principle for the soft terms, 
if one makes the assumption that gravity is flavor blind. SUSY is presumed to be broken 
at a scale Mj ~ 10 10 GeV in a "hidden" sector of particles which have only gravitational 
interactions with the "visible" sector of particles familiar to us. The SUSY breaking is 
transmitted to the visible sector by gravitational effects, so that the soft SUSY-breaking 
terms in the visible sector are characterized by a mass scale mguSY = -^j A^Planck ~ 
M% and are universal, independent of the unknown physics of the hidden sector. This 
mechanism implies a simple form for the soft breaking Lagrangian at the Planck scale. 
The following terms arise: 

• A common (mass) 2 (denoted TOq) for all scalars in the theory. 

• Scalar trilinear couplings which are given by the corresponding Yukawa couplings in the 
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superpotential multiplied by a common mass parameter Aq. 

• Scalar (mass) 2 terms given by the corresponding mass in the superpotential multiplied 
by a common mass parameter Bq. 

• A common mass m^ 2 for the gauginos. 

The mass parameters mo, Aq, Bq, and m^ 2 are all of order mguSY- The assumption 
of universality may have to be modified in non-minimal supergravity models [3] or in 
superstring models [4]; it is clear in any case that while they are sufficient to avoid flavor- 
changing neutral currents, they are certainly not necessary. The origins and consequences 
of scalar mass non-universalities which are safe for flavor-changing neutral currents have 
been discussed for example in [5-16]. 

In the MSSM, the superpotential is given by 

W = uY u QH u + dY d QH d + eY e LH d + (iH u H d 

where Q, L (it, d, e) are chiral superfields for the SU{2)i doublet (singlet) quarks and 
leptons; H U: H d are the Higgs doublet chiral superfields with weak hypercharge +1/2, —1/2; 
and Y u , Y d , Y e are 3x3 Yukawa matrices. The simple supergravity assumptions above tell 
us that the soft SUSY-breaking Lagrangian at Mpi anc k is 



L soh =- m l 



H u \ 2 + \H d \ 2 + (|g|2 + H 2 + | rf |2 + |L| 2 + |e|2) 

families 

Aq(uY u QH u + dY d QH d + eY e LH d ) + H.c. 



(1.1) 

- B fiH u H d + H.c. 

" \rn 1/2 {gg + WW + BB) + H.c. 

(Here we use the same symbol for the scalar fields as for the chiral superfields, and g, W, B 
are the gauginos for SU(3)c, 311(2)^, U(l)y respectively.) This parametrization serves as 
a set of boundary conditions on the model. The couplings so specified must be run down 
to ordinary energies using the renormalization group (RG) equations. This has been done 
by many groups; recent examples include [17-24]. 
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Of course, it is quite unlikely that there is no new physics between the TeV scale 
and the Planck scale. A much-heralded hint of this is the apparent unification of gauge 
couplings at Mjj » 2 x lO 1 ^ GeV with SUSY thresholds near the electroweak scale[25]. 
This suggests a SUSY grand unified theory (GUT) or superstring model, both of which 
in turn imply a variety of new chiral superfields and new gauge interactions at least two 
orders of magnitude below the Planck scale at which the boundary conditions (1.1) should 
be applied. The existence of this physics should be taken into account by running the gauge 
couplings of the theory from Mpi anc k down to the weak scale, integrating out non-MSSM 
fields along the way as they acquire masses. Unfortunately, it is very difficult to anticipate 
with any confidence the nature of the physics between Mjj and M-p\ &n( ±. For this reason, 
most studies have used the approximation of applying the boundary conditions (1.1) at 
Mjj rather than Mpi anc ^. Some recent papers[10,12,14] have taken into account the effects 
of RG running between Mjj and Mp^^, which can be large e.g. in certain GUT models. 

In this paper we will be mainly concerned with the D-term contributions [5] to soft 
scalar masses which should arise in any model with a gauge group of rank > 4 which 
breaks down to SU(3)c x SU(2)^ x U(l)y- Some of the phenomenological implications of 
D-term contributions have recently been explored in [6,7,9,10,11,15]. We will consider here 
the case that the underlying gauge group is an arbitrary rank 5 or 6 subgroup of Eq. This 
encompasses a very wide range of special cases, including SO (10) SUSY GUTs and string- 
inspired models including e.g. "flipped" 5*7(5) x U(l) and SU(3) C x SU(3) L x SU(3) R . 
The breaking of the additional U(l) factors embedded within such groups leads to specific 
patterns in the soft scalar masses which may be thought of as non-universal corrections 
to the universal boundary conditions (1.1), even though they typically arise at a much 
lower scale. These non-universal corrections to the scalar masses do not imply additional 
contributions to flavor-changing neutral current processes, since the contributions are the 
same for each set of squarks and sleptons with the same SU(3)c x SU(2)^ x U(l)y 
quantum numbers, e.g. (d^s^bi). The D-terms allow one to reach certain regions of 
parameter space which cannot be accessed using the universal boundary conditions, some- 
times with interesting implications for phenomenology below the TeV scale. For example, 
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much smaller values of \i can be obtained for fixed values of the other parameters if the 
D-term contributions have the appropriate sign and magnitude. The Higgs pseudoscalar 
mass might be lower than predicted in the universal case, and the lightest neutral Higgs 
eigenstate may have couplings which differ significantly from those of a standard model 
Higgs scalar of the same mass. The D-terms may also leave a dramatic imprint on sum 
rule relations among the squark and slepton masses. Such effects may eventually help to 
distinguish between various scenarios for gauge symmetry breaking at very high energies. 

This paper is organized as follows. In section 2 we will review the origin of D-term 
contributions to soft scalar masses, and introduce a parametrization of the most general 
such contributions arising in the MSSM from arbitrary subgroups of Eq. Using the RG 
equations we will show how the D-term contributions at very high energies leave their 
imprint on TeV scale physics. In section 3 we will discuss some of the salient phenomeno- 
logical consequences of the D-term-induced non-universality, taking into account direct 
and indirect experimental constraints. In section 4 we will make some concluding remarks. 



2. D-term contributions to soft scalar masses 

In general, D-term contributions to the scalar masses will arise whenever a gauge 
symmetry is spontaneously broken with reduction of rank [5] . To understand this, consider 
a toy model with a gauge group containing a product of abelian factors YljU(l)j. We 
take the U(l)j gauge couplings and charges to be normalized so that the U(l)j gauge 
supermultiplets have canonical kinetic terms. This determines, up to orthogonal rotations 
on the index /, a preferred basis for the U(l)j charges which is not invariant under general 
invertible linear transformations. The gauge group is to be spontaneously broken by VEVs 
for the scalar components of chiral superfields $ and $ with charges and — Q/<j> 
respectively under U(l)j. For illustrative purposes, we may for example suppose that the 
spontaneous symmetry breaking is accomplished by taking the superpotential to be 

W = \ — (n>2) 

with M a mass parameter of order the Planck or string scale. The supersymmetric part 
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of the scalar potential including D-terms is given by 

1 

" M 4n ~ 6 



vsusy =i(i $ P + m 2 )m 2n - 2 



+Ey(^(i $ i 2 -i $ i 2 )+E^i^i 2 ) 2 21 

I a 

where the additional fields (p a play the role of the scalar fields of the MSSM. In addition 
there are soft SUSY-breaking terms which include 

V soit = m 2 |$| 2 + m 2 |$| 2 . (2.2) 

The full scalar potential has a non-trivial VEV in a nearly D-flat direction: 

-(m 2 + m 2 )M 4 -- 6 l 1/(2n - 2) 



(2.3) 



(4w - 2) 

provided that m 2 + m 2 < at the scale (<&). This can be easily achieved if for example $ 
has large Yukawa couplings to some other superfields, driving the running m 2 negative in 
the infrared. The deviation from D-flatness is given by 

<d>> 2 - <$> 2 « \{m 2 - m 2 )/J2g 2 lQh • (2-4) 

After integrating out the fields $ and $, one obtains corrections to the soft scalar masses 
of the remaining fields <p a : 

AmJ = Qladi ; dj = \{m 2 - m 2 )g 2 jQ I(P / ^ 9jQ% ■ (2.5) 
I J 

From (2.5) we see that the D-term contributions do not depend on the precise form of the 
superpotential (2.1), and only depend on the soft masses m 2 , m 2 and the charges under the 
U(l)j gauge groups which participate in the breaking. The VEVs (2.3) reduce the rank 
of the gauge group by 1, with the surviving N — 1 i7(l)'s having charges which are linear 
combinations of the iV original U(l)j charges. The contributions to the (mass) 2 of the 
surviving scalar fields are just proportional to their charges under the original U(l)j gauge 
groups, with proportionality constants dj, and are always of roughly the same order as the 
original soft scalar (mass) 2 (i.e. to suSy)' even though the scale of spontaneous symmetry 
breaking set by (2.3) is many orders of magnitude larger [5]. 
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In a general situation, the rank reduction process exemplified above may be iterated 
several times; then it is clear that the total D-term contributions to the surviving scalar 
fields from the spontaneously broken U(l)j are of the form 

Am 2 a = J2QlaDi (2.6) 
/ 

where now the parameters Dj reflect the more complicated features of the symmetry 
breaking process, and in general parameterize our ignorance of such details as which fields 
are getting VEVs and possible hierarchies in the VEVs. 

In addition, if the spontaneously broken gauge groups included non-abelian generators 
T a which commute with all of the generators of SU(3)q x SU{2)j j x U{1)y, these could 
contribute residual soft mass terms of the form 

AL soft = - J D a (^t T » + H.c.. (2.7) 

However, in the MSSM we may dismiss this possibility if the underlying gauge group is 
a subgroup of Eq, or more generally if it does not mix families. This is because the only 
distinct chiral superfields with the same SU(3)c x SU(2)l x U{l)y quantum numbers in 
the MSSM are L and so that (2.7) would have to be of the form 

AL gofk = -D L L^H d + K.c. (2.8) 

which would necessarily imply R-parity violation. So we need only consider U(l)j D-terms. 
[It might be interesting to consider D-terms that could arise in variants of the MSSM with 
R-parity violation, or with family-dependent gauge symmetries (see e.g. [26]).] 

The group Eq contains two U(l) factors not contained in the standard model gauge 
group SU(S)c x SU(2)i x U(l)y- One may choose a basis for the two extra U(l) factors 
in an arbitrary way, but it is important to realize that the original U(l)j groups (with 
canonical kinetic terms) above the symmetry-breaking scale will in general be linear com- 
binations of the two chosen U(l) basis groups and the surviving weak hypercharge £7(1) y. 
We choose as a basis U(l)y, U(l)x, and U(l)g as in Table 1. 
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Table 1: U(l) charges of MSSM chiral superfields 





Q 


u 


d 


L 


e 


H u 


H d 


U(1) Y 


1/6 


-2/3 


1/3 


-1/2 


1 


1/2 


-1/2 


U(l)x 


-1/3 


-1/3 


1 


1 


-1/3 


2/3 


-2/3 


u(i) s 


-1/3 


-1/3 


-2/3 


-2/3 


-1/3 


2/3 


1 



One of the reasons behind this choice is that within a 27 of Eq there are two standard 
model singlet components v and S, which carry [£/(l)y, U(l)x, U(l)s\ charges [0, —5/3, 0] 
and [0, 0, —5/3] respectively. Thus U(l)x could be broken by a VEV for the scalar com- 
ponent of the field v (which carries lepton number —1), and U(l)g could be broken by 
the field S (which carries lepton number 0). The normalizations are chosen so that the 
largest charge appearing in Table 1 is unity; in order to obtain the correct normalization 
for unification into Eq (or any of its subgroups), these charges should be multiplied by 



3 3 3 
5' 2v / 10' 2VlO 



respectively. The gauged B — L subgroup of Eq is given by 

U(l) B -L = \u(l)y-hj(l) X . 



There is also a non-abelian SU(2) factor within Eq which commutes with all of the 
surviving SU (3)c x SU (2)^ x U(l)y gauge generators. It generates simultaneous rotations 
L <-> H^, v <-> S, and d <-> h, where h is a non-MSSM color anti-triplet. As we mentioned 
above, the off-diagonal generators of this SU(2) are prohibited by R-parity conservation 
from contributing to soft-scalar masses via D-terms. Its existence means, however, that 
there is a potential ambiguity in assigning U(l) charges to L, and d, stemming from 
different embeddings of these fields into remnants of a 27 for certain subgroups of Eq. The 
L <-> Hrf ambiguity is easily fixed just by the choice of U(l)x versus ^(1)5. [Note that 
under U(l)x <-> U(l)g, the charges of L and are exchanged and those of (Q, u, e, H u ) 
remain invariant.] The U (1) charges of d are then also uniquely fixed if we demand that that 
the MSSM Yukawa couplings are allowed by the unbroken gauge group. Other assignments 
of the U(l)x s charges for d would imply that B — L cannot be a gauged U(l) symmetry, 
a not particularly attractive possibility which we decline to pursue here. 



The D-term contributions to the soft scalar masses of the MSSM arising from sponta- 
neous breakdown of any subgroup of Eq may therefore be parameterized by: 

Am l=l D r-l D x-l D * a 2 jj 1 l D 

2 1 1 Am e = D Y ~ o D X ~ ~ D S; 

Am l = -~D Y - ~D X - -D S - ,12 2 

2 1 3 3 2 Am l, = + o D X + *D S ; (2.9) 
Ami = «D Y + D x - -De 



»j =3^+^-3^, 2 1 



Am| = - -£V + D x - -D S ; 2 6 

These corrections to the scalar masses should be applied at the scale Mjj of spontaneous 
symmetry breaking. In order to see how these contributions affect physics at ordinary 
energies, we must use the RG equations. For the squarks and sleptons of the third family 
and the Higgs scalars, these are given by 

dni ^ fi o 

16tt 2 ^ = -6^ 2 2 M| - -g\M{ - -g 2 S + 2y 2 E 2 T 

2 dm~ 24 2 2 , 6 2 C 1 /i 2^2 
d ^ 

16tt 2 ^^ = -6g%M 2 - jj^M? + ^g 2 S + 6y 2 E 2 

d'xxx^' (~\ ^ 

1 ^~£ ± = - 6 ^ M 2 " 5^1 M 1 " 5^ + 6 ^ E & + 2 ^ E ? 

where t = ln((5/<5o); -^3,2,1 are the running gaugino masses; yt,6,r are the running Yukawa 

couplings of the third family; A t ^^ T are the corresponding scalar trilinear couplings; S 2 = 

171% + m 2 + m 2 ~ +A 2 ;T l l = m 2 rr +m 2 + m? + Ah T, 2 = m 2 rr +m? + m\ +A 2 ;and 

S = Tr(Ym 2 ) = m 2 Hu - m 2 Hd + (m| - 2m| + m| - m| + m|) . 

families 

The gauge couplings (71 and cti are taken to be in a GUT normalization throughout 

this paper. The RG equations of the first and second family squarks and sleptons are the 
same except that the terms proportional to Yukawa couplings do not appear. 
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Now the quantities E|, E^ and E^ are unaffected when (2.9) are applied, and this 
persists as the scalar masses evolve according to the RG equations. So the only changes in 
the running of the soft scalar masses occur because of the presence of the terms proportional 
to S. One finds 

dS 66 a.\ 
~dt = ~~5 4^ 



which has the solution 



S(t) = S(t )ai(t)/ai(t ) . 



At the scale Mp associated with spontaneous symmetry breaking, we have from (2.9) 

S(t D ) = llDy + ^D X - l -D s 

Then at any scale lower than tjj we find that the change induced in the soft scalar masses 
by virtue of (2.9) (compared to a template model in which the D-terms are not present 
and all other parameters are the same) is given by 

Am 2 a (t) = Y a 



4(1 - r) „ (1-r) ' 



+ Qx«D x + QsaD S (2.10) 



with r(t) = cxi(t)/ai(tjj) and Y a , Qxai Qsa the charges of the scalar <p a under L r (l)y, 
U(l)x, U(l)g respectively. In general r decreases as one moves to lower scales; for 
near the GUT scale, r .43 near the electroweak scale. So, to a good approximation, the 
effect of RG running of the scalar masses is simply to reduce the Dy contribution by a 
factor of r, while leaving the Dx and Dg contributions untouched. If one is treating Dy, 
Dx, and Dg as parameters of our ignorance regarding the rank reduction mechanism, one 
may as well impose them at e.g. the GUT scale or at the electroweak scale. The error 
incurred in doing so can always be simply absorbed into a redefinition of Dy. We will 
henceforth refer to the effective value of Dy at the electroweak scale as Dy, it is equal to 
the quantity in square brackets in (2.10). 

In some special cases, it is possible to make statements about the relative sizes of Dy, 
D x , D s . For example, in an 50(10) SUSY GUT, in the limit that 50(10) breaks to 
SU(3)c x SU(2) L x U(l)y immediately at the GUT scale, one has 

D Y = D S = 0; D x ^0. 
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In the case that £0(10) breaks down to a smaller group [e.g. SU{A)q x SU(2) l x U(1)r] 
and only then to the MSSM far below the GUT scale, one has 



D s = 0; D Y ,D x ^0. 

However, there is another possible embedding of the MSSM into a SU{A)c x SU(2)i x U(l) 
subgroup of Eq which does not fit within £0(10); this would yield instead 

D x = 0; D Y ,D S ^0. 

In the case that the gauge group is a rank 5 subgroup of Eq obtained from a Calabi-Yau 
superstring compactification [27], one has 

D x = D S £ 

(with Dy possibly non-zero) because of the mandatory [28] presence of the subgroup 
U(l)x' = U(l)x + U(l)s in such models. It is useful to note that in general the magnitude 
of Dy is expected to vanish in the limit that U(l)y has a vanishing mixing angle with the 
spontaneously broken U(l) subgroups, a limit which is often approximately realized when 
the scale of spontaneous symmetry breaking is not far below the GUT scale. 

Experimental limits on sfermion masses imply constraints on Dy, Dx, and Dg. The 
condition that the right-handed selectron mass be above the present experimental limit 
amounts to approximately 

ml > (45 GeV) 2 (l - | cos2/3|) - .15m 2 /2 - Dy + -D x + ~D S , (2.11) 

while the lower limit on the mass of the sneutrino from its contribution to the invisible 
width of the Z boson is approximately given by 

ml > (41 GeV) 2 (l + 2.5| cos2/3|) - .5m 2 /2 + ^D Y - D x + ^D s . (2.12) 

(We assume tan/3 > 1 throughout this paper, so cos 2/3 < 0.) Either (2.11) or (2.12) 
typically sets the lower limit on m 2 , for given values of Dy, D x , Dg. 
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Note that if Dy — ^D x — g-Ds* and — ^D Y + D x — are both large and positive, 
one may entertain the interesting and unusual possibility that is significantly negative. 
Negative running squared masses for squarks and sleptons at high scales do not necessarily 
imply that the squark and slepton fields acquire VEVs, as long as there are no exactly 
D-flat and F-flat directions involving only MSSM fields in the supersymmetric part of the 
scalar potential at high energies. The correct place to study the scalar potential is at the 
scale of the putative VEV, which is typically much smaller than Mjj or Mpi anc ^; at this 
scale contributions to the running scalar (mass) 2 terms from gaugino loops will push them 
positive even if they were negative at very high scales, eliminating the purported VEV. The 
MSSM does have many directions which are exactly D-flat and F-flat at the renormalizable 
level. However, these can be lifted by non-renormalizable terms in the superpotential and 
by D-terms from additional gauge interactions which remain unbroken in the high energy 
theory. If for example U(l)g survives down to an intermediate scale Mj ~ ^/MzMjj [as 
suggested by (2.3) with n = 2], then all D-flat directions above Mj would have to involve 
VEVs for H u or H^, since all other chiral superfields in the MSSM have ^(l)^ charges 
of the same (negative) sign. One can show that all such flat directions can be lifted 
in models with e.g. chiral superfields and superpotential interactions for gauge-singlet 
neutrinos. Therefore, to maintain maximal generality, positivity constraints on linear 
combinations of scalar (mass) 2 terms which correspond to D-flat and F-flat directions of 
the MSSM at the renormalizable level should be imposed only at Mj, which we take to be 
10 10 GeV. Such constraints are easily incorporated into computer programs which calculate 
the sparticle spectrum numerically, as described in the next section. 

Even when the D-terms vanish, there is in the MSSM a possibility for m 2 , at the 
unification scale (or at least its effective value as seen from low energies) to be very slightly 
negative, but only if is large enough. Numerically, we find that in models which 

satisfy all other constraints, m 2 , > for m]/ 2 < 100 GeV, and ttiq > —(50 GeV) 2 when 
m l/2 < 400 GeV if D-terms are negligible. Also note that a negative value of Dg pushes all 
squark and slepton masses higher while pushing the Higgs (mass) 2 parameters lower. For 
this reason, we find that Dg has essentially no lower bound from experimental constraints 
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on sfermion masses alone, but turns out to be constrained by Higgs mass bounds and the 
stability of the Higgs scalar potential. 

Each of the contributions Dy, Dx, Dg will split m 2 H and 'm% d (in addition to the 
splitting induced by the top Yukawa coupling) so that there is a potentially significant 
impact on electroweak symmetry breaking. In particular, the value of \x needed for correct 
symmetry breaking (for other parameters fixed) may be raised or lowered. To understand 
this, we may consider the tree-level relation 

u 2 = -. — ^—^{m 2 Tr cos 2 3 - m 2 rr sin 2 8) - -Ml . 

P |cos2/3| v Hd Hu ' 2 z 

The change in /i 2 induced by the D-terms is therefore approximately 



A(^ 2 ) ' 



- l -D Y - 2 -D x + (cos 2 ii- 2 - sin 2 (3)D S 



(2.13) 



cos 2/3 1 

In models with universal boundary conditions and no D-terms, tends to scale with 
max(mo, mj m)- If the quantity in square brackets in (2.13) is negative, \fi\ can be sub- 
stantially smaller than is otherwise allowed, with important implications for the neutralino 
and chargino masses and mixing angles and for the Higgs scalar sector. 

The D-terms can also have a substantial impact on mass relations for the squarks and 
sleptons. For example, a sum rule given in [29] for the first and second family squarks and 
sleptons is modified to 

2m L - m \ R - m \ + m k - m L + y sin2 °wM 2 z \ cos 2/?| = -y D Y (2.14) 

Note that Dx and Dg only affect this sum rule indirectly (and weakly) through Dy, while 
a large contribution to Dy can cause a substantial deviation. 

Another interesting sum rule is [29] 

m\ R - m\ R = (.01)M| - | cos2/3|(43 GeV) 2 - D Y + ±D S - ^D x • (2.15) 

The first two terms on the RHS of (2.15) are generally small and of opposite sign, so 
that rriu R and mj^ are always very close to equal in models without D-term or other 
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non- universality; this relationship can be modified in a dramatic way by the D-term con- 
tributions as we will see in the next section. Using this and similar tests, a substantial 
deviation from the Dy = Dx = Dg = condition might be discernible at an e + e~ collider 
[30-32], yielding a strong (if not unambiguous) clue to physics at very high energies. 

Of course, D-term contributions will generally not be the only non-universal imprints 
left on the soft scalar masses by physics at very high energies. The very presence of 
additional gauge groups at high energies means that there will be additional contributions 
to the RG running. In addition, there will be the effects of running the RG equations 
between Mjj and Mpi anc ^. In order to keep our exploration of the parameter space within 
reasonable limits, we will neglect such effects insofar as they cannot be absorbed into the 
ignorance parameters Dy, Dx and Dg. This presumes that the additional gauge symmetry 
is broken not far below Mjj and that all interactions above Mjj are reasonably weak so that 
threshold and running effects are not overwhelming. This situation is preferred anyway if 
the apparent unification of gauge couplings is assumed to be not accidental. 

It may be useful, however, to introduce a parameterization of the most general family- 
independent scalar mass non-universality at Mjj. There is of course no unique way to pick 
such a parametrization, but we clearly want to arrange that the parameters Dy, Dx, Dg 
appear also in the more general parameterization, that RG running can be formulated in 
as transparent a way as possible, and that the parameterization be invertible in terms of 
the running scalar squared masses at the scale where they are presumed to be family- 
independent. Since there are 7 independent masses m~, m~, m?, m~, mjj , m|^), 
there should be 6 parameters (in addition to the common m^) which describe the non- 
universality at the input scale. We find it convenient to choose the contributions to scalar 
masses from the extra three parameters to be just proportional to those from the SU(3)c, 
SU(2)i, U{l)y gaugino loops to the RG equations for the scalar masses; we will call these 
contributions K$, K<i and K\ respectively. Thus one may parameterize the most general 
family-independent but non-universal scalar masses at the input scale (e.g. Mjj or Mpi anc k 
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in a given theory) as follows: 
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This matrix is invertible, so that in principle, if one were given the scalar (mass) 2 param- 
eters at the input scale, one could recover tuq, Dy^x,S an d -^3,2,1 : 
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Of course, in practice it is difficult to imagine being able to reconstruct the soft masses 
at the input scale, especially m? H and mjj . More importantly, the fact that the param- 
eterization is invertible means that it is also general. This parameterization has several 
other nice features. First, the i*Q contributions do not affect S, so that the running of the 
sfermion masses of the first two families are unaffected. Thus the imprint of non-zero 
on the squark and slepton masses of the first two families is exactly the same at the TeV 
scale as at the input scale. In particular, the sum rules (2.14) and (2.15) are not affected 
by non-zero K^s. At the low scale, the scalar masses of the first two families are given 
by an equation of the same form as (2.16), with Dx and Dg unchanged, and with the 
replacements Dy — > Dy and — > where Dy is just as before and 

K^Ki + Ciit); 



/3/5\ 
Ci(t) = 3/4 



tu 



dtaii^Miit) 2 , 



(2.17) 



\4/3/ n Jt 

taking into account one-loop RG contributions and neglecting the small Yukawa couplings. 
Non-zero K{S only enter the RG equations for the third family squarks and sleptons and the 
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Higgs scalars through the combinations T,j, and S^. These give additional contributions 
which may be described to a good approximation by the parameterization 

Am 2 iL ~ bL = -X t -X b ; Am 2 iR = -2X t ; Am^ = -2X b ; Am\ L ^ = -X T - 

Am 2 R = -2X T ; Am H u = - 3X f, Am H d = ~ 3X b ~ x t ■ 
When tan j3 is not large, X b and X T are negligible. 

In the MSSM, there is a large uncertainty in the numerical contribution of gluino loops 
to squark masses from the RG running, because the strong coupling constant is not very 
precisely known and because the function C^(t) in (2.17) runs quickly with scale below 
1 TeV. Note that this uncertainty can effectively be absorbed into the non-universality 
parameter K%. Similarly, the additional RG effects from MSSM gaugino loops (but not 
heavy non-MSSM gaugino loops) between Mjj and Mpi anc k, or above any intermediate 
scale in which the gauge couplings run differently (due to e.g. thresholds from heavy 
vector-like chiral superfields) can always be absorbed into the K^. Other new physics 
at very high energies (e.g. gaugino loops for gauge generators which do not commute with 
SU(3)(j x SU(2) L x U(1)y) may influence low energy physics through more complicated 
combinations of the parameters given above. 

3. Phenomenological constraints and numerical results 

In this section we will study some of the implications of D-term non-universality on 
SUSY phenomenology. To do so consistently requires that we simultaneously impose a va- 
riety of constraints from correct electroweak symmetry breaking, absence of color-breaking 
and charge-breaking global minima of the scalar potential, direct and indirect sparticle 
and Higgs boson mass limits [33], and a neutralino LSP. It is convenient to study the 
resulting constrained parameter space using a computer program which generates models 
randomly, imposing the constraints and calculating the resulting sparticle masses, mixing 
angles and couplings numerically. Here we describe the results of such an investigation, 
using a computer program similar to the one described in [22]. We first treat some spe- 
cial cases obtained by adding D-term contributions to a "template" model with universal 
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boundary conditions, to illustrate how experimental constraints limit the size of the D- 
terms and how the SUSY phenomenology is in turn modified by the D-terms. We will then 
turn to a much more general study of the parameter space allowed by arbitrary values of 
all parameters as constrained by experiment. In all cases, the D-terms are applied at a 
scale Mjj « 2 x 10 16 where the SU(3)c x SU(2) L x U(1) Y gauge couplings unify, and we 
retain the universal boundary conditions on gaugino masses and scalar trilinear couplings 
at Mjj. 

We choose as a template model the MSSM with the following parameters fixed: 
m\ = (100 GeV) 2 ; m 1/2 = A = 100 GeV; 

mtop = 175 GeV; tan/? = 3; sign( / u) = +; 

D Y = D X = D S = . 

Now consider turning on the D-terms. As explained in section 2, the following (approxi- 
mate) directions in parameter space may be theoretically preferred: 

(A) D x ^ 0; D Y = D S = 0, 

(B) D s ^ 0; D Y = D X = 0, 

(C) D S = D X ^0; D Y = 0. 

Therefore, we will investigate as examples what happens in each of cases A, B, and C 
separately. 

Case A: In Figure 1 we plot mg L , mg R , and as a function of D x , with D Y = Dg = 
and all other parameters fixed as in the template model. (In our notation, h and A are 
the lightest Higgs scalar and pseudoscalar boson, respectively.) The lower bound on D x 
in this example is set by the sneutrino contribution to the invisible width of the Z boson, 
which translates to a limit m e ~ L > 82 GeV in this example because of the sum rule 

m\ L -ml = |cos2/3|M^. (3.1) 

The upper bound on D x is set by the experimental bound [33] > 44 GeV. In the 
remaining range —(95 GeV) 2 < D x < (160 GeV) 2 , all other sparticles have experimentally 
allowed masses. Throughout most of this allowed range, is sufficiently large compared 
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to m/j so that h behaves essentially like a Standard Model Higgs boson. However, as 
Dx approaches its upper bound, and become small together, as does fx. This has 
potentially important implications for Higgs production via e + e~ — > Zh, the rate for which 
is proportional to the quantity sin 2 (/3 — a). [In the limit ^> m^, one has sin 2 (/3 — a) ~ 1; 
the heavier Higgs scalar eigenstates decouple.] For D x very close to its upper bound, the 
lightest Higgs mass eigenstate h has couplings which are not like a Standard Model Higgs 
boson; for Dx > (156 GeV) 2 , we find in this example sin 2 (/3 — a) < 0.6. However, as long 
as D x < (135 GeV) 2 , one has sin 2 (/3 - a) > 0.9 here. 

Case B: In Figure 2 we show mg L , mg fl , and this time as a function of Dg added to 
the template model. As can be clearly seen from (2.9), mg L , m^ R and all other squark and 
slepton masses grow monotonically with increasingly negative Dg. Thus the lower bound 
on Dg is typically set by the lightest Higgs mass bound; here we find Dg > —(325 GeV) 2 . 
Conversely, the upper bound on Dg is usually set by the lower experimental bound on either 
mg fi or (as in this case) m^. For Dg near its experimental lower bound, h behaves very 
differently from a Standard Model Higgs, with sin 2 (/3 - a) < 0.6 for Dg < -(310 GeV) 2 . 
However, as long as D5 > —(260 GeV) 2 , we find sin 2 (/3 — a) > 0.9 in this example. 

Case C: In Figure 3, we again plot mg L , m^, and m^, now as a function of Dx = Dg 
with Dy = and all other parameters fixed as in the template model. In this case, nothing 
very remarkable happens to the Higgs masses or couplings over the entire allowed range for 
Dx = Dg. Both mj{ and fall somewhat as Dx = Dg increases, but sin 2 (,9 — a) > 0.96 
always in this example. The lower bound on Dx = Dg is set by the sneutrino contribution 
to the invisible width of the Z boson (or, equivalently, by mg L > 82 GeV for this value of 
tan/3). The upper bound on Dx = Dg is set by the lower limit of 45 GeV on mg R . Note 
that for Dx = Dg less than only —(59 GeV) 2 , one finds m^ R > mg L , in contradistinction 
to the usual situation in models with universal boundary conditions. 

The rate for Higgs scalar boson production from e + e~ — > Zh is equal to sin 2 (/3 — a) 
times the corresponding rate for a Standard Model Higgs scalar boson of the same mass. 
The situation (illustrated by cases A and B) of light pseudoscalar masses together with 
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small values of sin 2 (/3 — a) can occur for all values of tan/3 and mtop when appropriate 
D-terms are present, unlike in the case of universal boundary conditions. For all values of 
tan/3, we found some models for which sin 2 (/3 — a) was very close to zero. This generally 
happens near the "edge" of D-term parameter space, i.e., for values of Dy, D x , and Dg 
near their limits. When this occurs, the production cross-section for h is reduced, but 
the rate for hA production is correspondingly enhanced, so that hA production may be 
observable at LEP2. 

We now consider the possibility of arbitrary values for all of the parameters. The re- 
maining results of this section were obtained by an exhaustive exploration of the parameter 
space, with the following constraints: 

-(200 GeV) 2 <ml< (800 GeV) 2 ; 
40 GeV < m 1/2 < 400 GeV; 
1.5 < tan/3 < 60; sign(^) = ±; (3.2) 

160 GeV < m t0 p < 190 GeV; 
-(500 GeV) 2 < D Y ,D X , D s < (500 GeV) 2 . 

As explained in the previous section, we allow for the possibility of negative m 2 ). (We 
have checked, however, that this possibility does not have a strong impact on the general 
results outlined below, although it clearly can affect physics within individual models.) We 
impose constraints on negative squark and slepton running (mass) 2 terms at Mj = 10^ 
GeV corresponding to each of the flat directions (w, e); (w, <i); (Q 7 L); (<i, L); and (e, L) (in 
each case with appropriate flavor structure). Here (u,e) denotes for example the D-flat 
and F-flat direction 






V2v , 



which implies the constraint 2m? + m? + 2m| > at Mr. 

^ Ur t R T R 1 



The range allowed for Aq was determined for each model by the requirement that there 
be no color-breaking global minima of the scalar potential at the electroweak scale in the 
D-flat (but not F-flat) directions (t L ) = (tR) = (H u ) ^ and (b L ) = (b R ) = (H d ) ^ 0. 
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In addition, we required that there not be global color-breaking minima in the non-D-flat 
directions (ti) ^ or (£#) 7^ 0, in each case with all other VEVs vanishing. The latter 
requirements amount to approximately 

mj > -Ml\cosm\( a3 + 3a2/4 ) 1/2 (3-3) 
*£ Zl V n V3a 2 + 9ai/5y V ' 

m? > -Ml|cos(2/3)|f Q3 + 4ai/ f (3.4) 

respectively at the weak scale. Here m~ and m~ are the running stop mass parameters 
in the absence of electroweak breaking, which appear in the usual stop mass matrix as 

m i L + m top ~~ - 35 l cos 2/5 1 M§ mto P (A t + [i cot (3) 

mt ov {At + // cot/3) m?^ + mj 0p - .15| cos 2/3 |M| 

The constraints (3.3) and (3.4) effectively eliminate the possibility of a stop squark much 
lighter than mtop with a small stop mixing angle, and thus limit certain combinations of the 
D-terms. In this study we do not allow the Ki introduced in the previous section to vary 
because of the practical need to confine the scope of the investigations, although it would 
certainly be interesting to consider this. We should remark that the limits on the dimen- 
sionful parameters in (3.2) were motivated both by practical considerations (e.g. computer 
running time) and our theoretical prejudices. Some of the results below would have to be 
generalized if, for example, we allowed jtiq < —(200 GeV) 2 , or larger D-terms. With that 
caveat, we now proceed to see how the allowed supersymmetric spectrum is modified by 
the presence of D-terms. 

Let us first consider the fate of some correlations between sparticle masses which may 
be taken for granted in the case of universal boundary conditions. 

In Figure 4, we show the allowed region in the (mg L ,mg ii ) plane for all models with 
universal boundary conditions (bounded by solid lines) , and the additional region (bounded 
by dashed lines) allowed for all models with arbitrary D-terms. In the universal case, viable 
models lie within a fairly restricted wedge; in particular, mg L > va^ R holds as a general rule. 
(If we relaxed the restriction TO]y 2 < 400 GeV in our original choice of parameter space, 
the region allowed by universal boundary conditions would be enlarged for mg L > 280 

20 




GeV.) With D-terms present, and mg L can each independently essentially saturate 
their lower limits. The lower limit on mg L in all models considered here is about 75 GeV; 
except for small tan/3, this limit comes from the sneutrino contribution to the invisible 
width of the Z and the sum rule (3.1). As we saw from our examples A, B, and C, the 
slepton masses often (but certainly not always) set the limits on the D-terms. 

Figure 5 depicts in the same way the allowed regions in the {m^ L ,m^ L ) plane. Here 
the correlation in the universal boundary condition case is looser, because m^ L grows 
much faster with m-j/2 than mg L does. In the presence of D-terms the allowed region is 
considerably larger, but not quite enough to saturate the experimental lower limits on 
squark masses when mg L > 150 GeV. Still, a substantial area is opened up for which 
m^ L < mg L , a situation which does not occur in the case of universal boundary conditions. 
A similar inversion of the usual mass hierarchy can occur for each of the other pairs of 
squarks and sleptons, depending on the values of Dy, Dx, and Dg. 

The masses of the first and second family squarks are always very strongly correlated 
with each other in the case of universal boundary conditions. In Figure 6, we show the 
(extremely narrow) allowed region in the (m^ fl , ) plane with solid lines, and the much 
larger region allowed by D-term contributions is enclosed in the dashed lines. However, 
the correlation between the two squark masses is not completely destroyed by D-terms, 
in the sense that e.g. m^ R < 200 GeV implies mj^ < 300 GeV; the converse is not true, 
however. The squark masses do not always saturate their lower bounds from direct searches 
(although mj^ can be quite small), mainly because the relative magnitude of the D-terms 
is constrained by the lower limits on the slepton masses, by Higgs boson searches, and by 
correct electroweak symmetry breaking. 

Figure 7 depicts the allowed regions in the (m^,m^) plane. The wedge (bounded 
by solid lines) allowed by universal boundary conditions is slightly wider than in Figure 6. 
This time the correlation between the squark masses is much stronger than in Figure 6, even 
for D-terms as large as allowed by other constraints. This can be easily understood from 
(2.9); since m\ R and w?~ obtain the same contribution from Dx and Dg, a big difference 
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between them can only occur if Dy dominates. However, such a situation is limited by the 

lower bound on mg R (if Dy > 0) or by the sneutrino contribution to the invisible width 

of the Z boson (if Dy < 0). There is an interesting region for which m^ R > mj , m^ L ; in 

this case the number of charginos and second-lightest neutralinos from gluino decays might 

be substantially increased over the usual expectation from universal boundary conditions. 

This can also occur if m,- > m,- , m,-,. . 

d R d L ' u l 

The mass of cIr (and b\) can even be much less than the lightest stop when D-terms 
are present. In Figure 8, we show the allowed regions in the (mj ,m^) plane. It should 
also be noted that with D-terms present, cIr especially can be considerably lighter than the 
gluino, a sharp contrast from the usual rule in the case of universal boundary conditions 
that all first and second family squark masses are greater than about .85 Mg. 

As remarked in section 2, large D-terms can also have a significant effect on the neu- 
tralino and chargino sectors, through their impact on the \i parameter as implied by (2.13). 
In much of the allowed parameter space in the MSSM with universal boundary conditions, 
li is large compared to so that the LSP (N\) : the second-lightest neutralino (iV^), and 
the lightest chargino {C\) are usually gaugino-like. In Figure 9, we show the allowed re- 
gions in the (p, M<i) plane for universal scalar masses (within solid lines) and for non-zero 
D-terms (within dashed lines). In the latter case there is a substantial new region with 
relatively smaller \i. The region with \i 3> M<i is also enlarged. If we restricted ourselves 
to any particular fixed values of tan/? and mt p, the enlargement of the allowed region 
in the (jU, M<i) plane would be even more dramatic. We find that eq. (2.13) is usually a 
quite good indicator of the effect of the D-terms on fx, even with full one-loop radiative 
corrections in the Higgs sector taken into account. 

In the MSSM, there is a well-known correlation between the masses of the lightest 
chargino and the second- lightest neutralino. In the limit that M^/(//±Mi 2) is small, one 
finds [29] 

rrir, ~ m Mo s ^ 

However, when ~ M2, this expansion fails so that one might expect a rather different 
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relation between these masses. In Figure 10 we show the allowed regions in the (m^v , ) 
plane for universal and for D-term boundary conditions. The latter region is significantly 
larger, particularly for smaller fj,, M^. On the other hand, we find that other important 
features of the chargino and neutralino sector are largely unchanged by the presence of 
D-terms. An important example concerns the mass differences ~ m N l an d m N 2 '~ m N l 
which are important in determining the visibility of a possible trilepton signal for SUSY 
at a hadron collider, or of chargino pair-production at an e + e~ collider. We find that 
~ m N 1 > 15 GeV always within our parameter space (3.2), even when is lowered 
with respect to Mi by D-terms. The allowed regions in the (m ^ , ) plane are shown 
in Fig. 11 for universal (enclosed in solid lines) and D-term (enclosed in dashed lines) 
boundary conditions. The charged lepton coming from chargino decay C\ — > N\F~v should 
thus always be energetic enough to serve as a trigger, if the corresponding branching ratio 
is large enough. The mass difference m ff 2 ~ m fi l ls also always larger than 7 GeV (and in 
the vast majority of models we found it was larger than 15 GeV). In Figure 12 we show 
the allowed regions in the ( m jf i m N ) plane. Of course, these results are really more of a 
statement about the form of the neutralino and chargino mass matrices (assuming gaugino 
mass unification) after all present experimental constraints have been imposed than it is 
about the consequences of D-terms, which only affect this sector indirectly through their 
effect on /x. 

[In the case that fx <C m>i/2-> one could imagine having almost degenerate C\, N21 
and N\ with large higgsino content. While this is a logical possibility, it is impossible 
to achieve in practice without extreme fine-tuning of parameters including Dy, R>x-> and 
Dg. In particular, it did not occur in the search of parameter space described above. This 
situation would also imply that the relic density of LSP's would be far too small for them 
to constitute the cold dark matter.] 

We conclude this section by mentioning the effect of D-terms on the Z-peak observable 
Rfj = F(Z — > bb)/F(Z — > hadrons). The Standard Model prediction for this quantity is 
lower by about 2a than the experimental value = 0.2202 ± .0020 obtained at LEP[34]. 
Supersymmetric corrections to are almost always tiny [35-38]. In particular, in the 
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MSSM with universal boundary conditions, Rfj never gets closer than 1.5a to the central 
experimental value. In the most general supersymmetric scenario, the prediction for 
can only be reconciled with the la experimental range if certain strong conditions are 
met. For tan/3 < 30, these include light stops which are predominantly £r, and light 
charginos with a significant higgsino content. One might imagine that these conditions 
could be attained in the presence of D-terms. Somewhat surprisingly, we find that this is 
not the case; an examination of viable models satisfying the constraints (3.2) finds that the 
prediction for can only slightly exceed the values indicated in Figure 1 of ref. [36] , and 
never enter the la experimentally allowed range. This may be understood qualitatively 
from the previous results; when the superpartners are near their experimental lower bounds 
(as required for larger R^), one finds that only relatively mild deviations from universal 
boundary conditions can occur. The extent to which the mass of the lightest stop squark 
can be lowered without having a large stop mixing angle is limited by the constraint (3.4) 
on color-breaking minima. On the other hand, a large stop mixing angle generally implies 
that is relatively large (so that the lighter chargino is not higgsino-like) or that At 
is large (which is problematic because of possible color breaking minima of the scalar 
potential in a D-flat direction). This implies that the experimental value of R^ must fall 
if minimal SUSY is correct with universal gaugino mass boundary conditions and only 
D-term contributions to scalar mass non-universality as imposed here. 



4. Conclusion 



In this paper, we have examined the possibility of D-term contributions to soft scalar 
masses in the MSSM which arise if the unbroken gauge group at very high energies is 
an arbitrary subgroup of Eq. These contributions are parameterized by three quantities 
Dy, Dx, and Dg. The effect of one- loop RG evolution on these parameters can always 
be absorbed into a redefinition of Dy. The sizes of Dy, Dj, and Dg (for a given set of 
other parameters) are limited by direct and indirect constraints on sparticles (especially 
sleptons) and Higgs scalar boson masses, as well as from correct symmetry breaking. 
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Of course, there is no guarantee that the D-terms will be large enough to make their 
presence felt in comparison to the universal contributions to scalar masses from mg and 
gaugino loops. Furthermore, it is clear that a future observation of apparent scalar mass 
non-universality at Mjj may not be unambiguously ascribed to D-terms, without further 
theoretical input. Still, checks of relationships which hold in the MSSM between squark 
and slepton masses will provide an important clue to physics at very high energies. With 
experimental constraints imposed, we found that the impact of the D-terms might be most 
dramatic on squark and slepton mass relations [see, for example, eq. (2.15) and Fig. 7]. 
In several important respects, we found that deviations from the universal boundary con- 
ditions tend to be rather mild, especially when superpartners are relatively light. For 
example, we did not succeed in finding a viable model with non-zero D-terms for which 
was within la of the present experimental value. We did find models at all values of tan (3 
and mtop for which sin 2 (,9 — a) was close to zero, in which case Higgs production through 
the usual Standard Model mechanism e + e~ — > Zh will be impossible to detect. How- 
ever, this generally occurred only when the D-terms were quite close to their experimental 
bounds, except for large tan/?. If scalar mass non-universality proves to be a phenomeno- 
logical necessity, D-term contributions will be an attractive theoretical candidate, since 
they maintain the natural suppression of flavor-changing neutral currents. 
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Figure Captions 

Figure 1: The masses of the lightest Higgs scalar boson h (lower solid line), pseudoscalar 
Higgs boson A (upper solid line), left-handed selectron (dashed line), and right-handed 
selectron (dot-dashed line), as a function of the parameter Dx added to the template 
model described in the text. 

Figure 2: The same quantities are plotted as in Figure 1, but now as a function of Dg 
added to the template model described in the text. 

Figure 3: The same quantities are plotted as in Figure 1, but now as a function of 
Dx = Dg added to the template model described in the text. 

Figure 4: The allowed regions in the (mg L ,mg fi ) plane for universal boundary conditions 
(inside solid lines) and arbitrary D-term boundary conditions (inside dashed lines), for 
models satisfying (3.2). 

Figure 5: The allowed regions in the (mg L , m^) plane. 
Figure 6: The allowed regions in the (m^, m J fl ) plane. 
Figure 7: The allowed regions in the (m^, m J L ) pl ane - 
Figure 8: The allowed regions in the (mj , rn^) plane. 
Figure 9: The allowed regions in the (/i, M2) plane. 
Figure 10: The allowed regions in the [m^ ,rn^) plane. 
Figure 11: The allowed regions in the (m^ , tUq ) plane. 
Figure 12: The allowed regions in the (m^ > m jy 2 ) plane. 
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